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Abstract 

We show that an electronic beam passing through the colhsionless plasma of the "cold" ions and 
the "hot" Boltzmann electrons can give rise to the propagation of the supersonic ion-acoustic rarefac- 
tion shock waves. These waves are analogous to those predicted by Zeldovich |5| in gasodynamics 
and complementary to the ion-acoustic compression shock waves in collisionless plasma described by 
Sagdeevia. 
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1 Introduction 

It is well known that the supersonic compression shock (CS) waves can propagate in the ideal gas, but 
their alternative - the supersonic rarefaction shock (RS) waves - cannot. In 1937 Zeldovich (51 has shown 
that the interaction between the gas molecules possessing the repulsive and attractive parts (e.g., the Van 
der Waals' gas) results in substantial changes: the only RS waves could been stably propagated, but the 
CS waves are unstable. 

In the 1960th Sagdeev O has predicted the CS wave which may propagate in the collisionless (CL) 
plasma. This appears when the CL plasma comprises the "cold" ions and the "hot" Boltzmann electrons 
and the mean free path is determined by the Debye radius of electrons. The shock wave is arisen by the 
supersonic compression soliton disturbed by dissipative effects. Such CS waves are analogous to those 
observed in gasodynamics. 

Regarding the RS waves in the CL plasma, a study of such waves dates back to the papers [1] and H 
for the media comprised the electrons with bi-Maxwellian distribution function and "cold" ions. 



In the present article we suggest another kind of the RS waves generation in the CL plasma based on 
the supersonic rarefaction solitons. We consider a problem in two steps. First, we show that the supersonic 
rarefaction soliton (of density and electric potential) can propagate in the CL plasma supplemented by 
the electronic beam. Next, being partially reflected by the soliton's electric potential, the beam disturbs 
its symmetric shape in such a way, that behind the maximal value of potential there appear nonlinear 
oscillations. 

2 Setup of the Problem 

Consider the CL plasma with the Boltzmann distribution Nf. = N^gexp {e(l){x) / nTf.) of the "hot" elec- 
trons, where N^o is a concentration of hot electrons in homogeneous plasma, Tg is a temperature of elec- 
trons, e and k stand for the charge of electron and the Boltzmann constant, respectively. The electric 
potential is varying along the x direction together with the charge density. Such plasma with elec- 
tronic beam can be realized in experiments with a hollow anode plasma source 13. In what follows, we 
use the renormalized potential il){x) = —\e\(t){x) / kT(. and choose the reference frame related to the elec- 
tric potential and normalized coordinate ^ = x/D where D = ^Je^^T^jNi^ denotes the Debay radius 
and Nio and eo stand for homogeneous ions density and dielectric vacuum permittivity, respectively. In 
this reference frame the ions are running with velocity Vi^) toward the wave front which is defined by 
requirement V when '0 — )• 0. Here U denotes the velocity of the homogeneous flow of plasma. 

Equations of the energy and mass conservations for cold ions describe their 1-dim stationary motion 
and lead 131 to the non homogeneous ions density Ni{^), 



where rrii and C denote the mass of the ion and the velocity of the ion sound. The Max number for the 
ions motion is denoted by M . 

Let £ = mv"^ / {2kT(,) + ip{C) be an energy of the single electron of the beam in the kT^ units. Then 
the density Njj{ip) of the beam electrons is given by the distribution function /(^, v). 



where Ni,o denotes a density of the electronic beam in homogeneous plasma, and m and v stand for mass 
and velocity of electron. 

The dimensionless potential Tp{^) is a smooth function and obeys the Poisson equation. 
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where 

BiC) = ^, B{0) = ^, A=§^, A + BiO) = l, A<1. (4) 

The last equahty holds due to the quasineutrality of homogeneous plasma. 

Regarding Ni,{il)), in this paper we deal with two different types of the density distribution of electronic 
beam in accordance with two posed problems: the existence of the supersonic rarefaction soli ton (see 
section |3]l and appearance of the supersonic RS waves (see section |4]l. 



3 The Rarefaction Soliton 

Consider the electronic beam with a distribution function f\E) given as follows, 

f{e) > if ^^8i<8 <£2 and f{e) = if ^: < or ^: > ^2 • (5) 
Then, by dU, @ and ^ we have B{S,) ~ ^(0), and equation ([3]) can been integrated as follows, 

^(^)' = M2(^y^7^-l^+^(e~^-l)-i?(0)^ + i^, K = const. (6) 

Equation Q gives rise to soliton solution if dip /d^ — when ^ — )• 0, that immediately requires K = 0. 
Rewrite equation ^ for tp <^ I and preserving the ip- and ip"^ terms, 

^(|)^,a-..«(o„4(.-i,). 

The linear in ijj term in ^ is vanishing due to and what is left results in M > 1/^/A. Keeping in mind 
A < 1 we get the supersonic rarefaction soliton, M > 1. 

Since 'ip{^) is a smooth function everywhere (also when ip > I) and ip{^) arrives its maximal value 
Tpmax then by Q we get 

^1 + ^^ -l^+A (e-'^-^ - l) - Bmma. = . (8) 

The last equation can be resolved analytically as M = M{A, B{0),'ipmax)- E.g., for the experimental date 
ia B{0) ~ 0.3, A ~ 0.7 and ip max — 20 we obtain M ~ 1.035. In Figure [T] we present the phase 
portrait of equation ([T]) for the parameters given above. The center (the singular point) is corresponded to 
K ~ 0.579 while the separatrix (soliton) is corresponded to K = 0. 

Up to the sign of the potential ip{£,) the phase portrait at Figure [T]is similar to that of the nonlinear ion- 
acoustic travelling waves of compression [3]. The main difference appears: in the latter case the maximal 
value of the ion density Ni{^) and the maximal value of the potential ip{£,) occur at the same coordinate ^, 
while in the former case the minimal value of the ion density occur when the potential arrives it maximal 
value. This is why the soliton solution ip{^) is accompanied by the ionic rarefaction. 
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Figure 1: Phase portrait of Eqn ^ for travelling waves of rarefaction; a soliton is drawn in bold. 

4 The Rarefaction Shock Waves 

In the model of the CS waves |]3l a symmetric shape of soliton is disturbed when a small portion of ions 
(with a lower energy) is reflected by potential barrier, but a large portion of ions (with a higher energy) is 
passed throughout it. In our case the rarefaction soliton will be disturbed by electrons of the beam which 
are reflected by the barrier. 

Choose the step-like distribution function f{£), 



m = r/'" e{£2 -£)&{£- £1) , (9) 

V uTe ^£2 - \J£\ 

where 6(a) denotes the Heviside step function. Substituting Q into Q we get a relative density of the 
transmitted electronic beam, 

Formula (flOl ) holds for the monotone growing potential -0(0 when ^/;(^) < £\ < £ < £2, ox ^ < ^1, 
where V'(Ci) = £\ (see Figure|2l). Here the incident electrons do not reflected by barrier but are transmitted 
if£: > £x. 

Next, let us focus on the other case, £\ < ^{(^ < £ < £2, where V'max = fpi^max) and ^max stands 
for location of maximal soliton potential. Here the beam is partly penetrated into the soliton potential 
when ^1 < ^ < Cmax, but the rest of electrons are reflected from the barrier when ■ip{£,) > £1, or ^ > ^1. 
It leads to the changes in formula (ITOl) 



M0 = By^£2-i'{0 , ei<e<W, B= J^^°^^ , (11) 
Finally, consider the last case i^max < £ < £2 (electrons do not reflect by the barrier) 



For convenience, unify three functions i?t(0, i = 1,2, 3, given in different ranges of ^ by one /(^) 
given in wiiole range of ^, 



- Cmax) Vi'max " V'(O) • 



(13) 



In the similar way we can construct the relative density R{0 of the electronic beam reflected from the 
barrier, 



R{0 = B [v^max - m - (^1 - m) V£i - i^io) , e < c« 



(14) 



A choice of arguments in the 0-functions in ([T3] ) and ([T4l ) is motivated by need to represent the Poisson 
equation Q as an autonomous differential equation. 
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Figure 2: The electric potential '0(C) for the nonlinear RS wave. Its 1st maximum ipmax satisfies inequality 

max <£2. 

Substitute the entire density B{^) = + R{C) of the beam into equation ^ and obtain its 1st 
integral, 

2 



1 fdip 

where the quasipotential V- {if)) is given by 

M2 



y_(^)=o, i<ir^ 
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(^:2 - V^)3/2 + {i^max - ^f'^ - 2(^1 - V)'/'0(^i - V') + 2ff/' - sT - 



In (fTSl ) we have taken a zero's value for the integration constant (see (O with K = 0) to provide the 
requirement dijj /d,^ — when ifj ^ {). 



Rewrite the quasipotential ([T6l ) for ^ <C 1 and preserve the non-linear terms up to ijp, 
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(17) 



Note that by (I10I14I) the term B ( + 



which enters into the linear in V' term in ([TT] ). gives a total density of the electronic beam including inci- 
dence and reflection as well. Then the whole linear in ^ term in ([TTl) disappears due to the quasineutrality 
of homogeneous plasma in the general case (when beam's reflection exists), 

(18) 



A^B (v^+ VVmax - 2v^) = 1 ^ A<\. 



Substituting the quadratic in -0 term of dlT] ) into ( [T5] ) we find the lower bound for M, 
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(19) 



Combining (fTSl ) and ( fT9l ) we arrive at M > 1/A > 1, i.e., the supersonic RS wave. 

An exact value of the Mach number can be found if we consider equation ([TSll for ^ = ^max and 

= V'max > "fi, i.e., when di)/di = 0, 
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(20) 



The last equality together with a quasineutrality condition (1181 ) allow to obtain M and A if the other four 
parameters B, £\, £2 and ipmax are given. 

Now, consider the Poisson equation Q in the range ^ > imax 

2 

- y+ (V') = , ^ > Uax , 



dip 



(21) 



where the quasipotential V"+(V') is given by 
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(22) 
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{£2 - iPf'^ - ii^max - ^f'^ - {£2 - i^maxf' 



In (I22I) we have taken non-zero's value for the integration constant (see ^ with / 0) to provide 
the requirement Vj^{'pmax) = 0. Note that if the potential ^(^) arrives also its minimal value Vmm at 
iniin > Cmax (scc FigurelJ]), then V"+(V'mm) = as it follows from (|2T]) . By equalities (ITSl) . (l20l ) and (l22l) 
it follows also that the value ipmin is completely determined by four parameters B,£i, £2 and ipmax- 
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Figure 3: Phase portrait of equations (1151) and (|2T)) for the RS waves. 

In Figure [3] we present the phase portraits of two first integrals (ITSl) and (|2TI ) of the Poisson equation 
Q for the date £i = 15, <?2 = 20, ipmax — 16.5 and B{0) ~ 0.15 taken from experiments 121. The values 
A ~ 0.803, M ~ 1.433 and iprnin ^ 12.86 were found by ^ and respectively. 

Both portraits are unified in one smooth curve. Its 1st part (an arch in the upper half -plane, dip/d^ > 0) 
is related to the soliton-like behavior of the electric potential ^p at the front of the RS waves, while its 2nd 
part (a closed loop) is related to the nonlinear oscillations of t/j. 
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Figure 4: Densities' distribution of the total electronic beam -B(^), including reflections (plain); of the 
transmitted electronic beam I(^) {dashed), and of the ions G{£) {bold). 

In Figure |4] we present three densities' distributions: the total density of the electronic beam B{£) 
when the reflected electrons are accounted for; the density of the transmitted electronic beam /(^) and 
also the density of the ions in plasma, G{£). The maximal value of G(^) is arrived before the wave front 
{^ = 0) of electric potential "0(0' while the density of the ions in the region ^ > is much less. This 
manifests that we get the rarefaction shock wave. 

We appreciate useful discussions with L.P.Pitaevsky, Y.E.Krasik and I.D.Kaganovich. 
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